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Abstract. Single-particle dynamics of the Anderson impurity model are studied 
using both the numerical renormalization group (NRG) method and the local moment 
approach (LMA). It is shown that a 'two-self-energy' description of dynamics inherent 
to the LMA, as well as a conventional 'single-self-energy' description, arise within 
NRG; each yielding correctly the same local single-particle spectrum. Explicit NRG 
results are obtained for the broken symmetry spectral constituents arising in a two- 
self-energy description, and the total spectrum. These are also compared to analytical 
results obtained from the LMA as implemented in practice. Very good agreement 
between the two is found, essentially on all relevant energy scales from the high-energy 
Hubbard satellites to the low-energy Kondo resonance. 



1. Introduction 

The numerical renormalization group (NRG) method [1, 2] forms a well established, 
powerful numerical technique for calculating the properties of quantum impurity and 
related models. The basic paradigm here is of course the celebrated Anderson impurity 
model (AIM) [3], in particular the Kondo effect arising in strong coupling (see [4] for 
a review). NRG is not restricted to calculation of static properties, but can equally 
handle dynamical properties such as single-particle spectra, see e.g. [5-7]. 

Dynamics in particular pose well known difficulties for analytical approaches [4]; 
be it in handling strong correlations in general, spanning the full relevant range of 
energy scales, satisfying the low-energy dictates of Fermi liquid behaviour and so on. In 
recent years we have been developing a local moment approach (LMA) to single-particle 
dynamics and related properties of a range of quantum impurity models, see e.g. [8-14]. 
Physically transparent, and technically straightforward in practice, the LMA is based 
on an underlying 'two-self-energy' description. As in Anderson's original work [3], local 
moments are introduced explicitly from the outset, leading to two degenerate mean- 
field saddle points. Spin-flip tunneling between the mean-field states — embodied in 
dynamical contributions to the two associated self-energies, naturally absent at pure 
mean-field level — lifts the erstwhile spin degeneracy of the saddle points, and restores 
the singlet symmetry characteristic of the local Fermi liquid state. 
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As implemented in practice the LMA is of course approximate, as dictated by 
the particular approximation chosen for the dynamical contributions to the two self- 
energies [8-13]. But it is taken for granted that, in principle, dynamics may be obtained 
either within an underlying two-self-energy framework, or within the 'single-self-energy' 
description that provides the conventional route to single-particle dynamics. In that 
case, it is natural to ask whether a two-self-energy description arises also within an 
essentially exact numerical approach such as NRG. It is this issue we consider here. 
We show that according to whether even or odd RG iterations are considered, both 
the single- and the two-self-energy descriptions are contained in the NRG approach; 
and in particular we obtain explicit NRG results for the composite broken symmetry 
spectra inherent to a TSE description. That in turn enables direct comparison to 
be made between NRG results for the underlying spectra and those arising from the 
LMA as implemented in practice. Very good agreement is found, on all energy scales 
characteristic of the problem. Results are given in §3, following a brief discussion of the 
background theory in §2. 

2. Theory 

The AIM Hamiltonian [3] is given in standard notation by 

H = ^ ek^ka + 5^(ei + ^ni--a)nia + Vikicj^c^^ + h.c.) (2.1) 

k,cr (T k,<T 

where the first term refers to the non- interacting host, and the second to the impurity 
with local interaction U and energy e^. In strong coupling (large U) the low-energy 
physics of the AIM is of course that of the Kondo model [4]. Our focus is the T = local 
impurity spectrum D{uj) = -ilmG'(u;), with G{uj) (^ G{t) = -iO{t){{Ci^{t),c\^})) the 
impurity propagator, 

G(uj) = T^- (2.2) 

Here g{oj) = [uj~^ — — A(a;)]~-^ is the non-interacting propagator with A(cj) = 
Sk^kl'^''' ~ ^k]~^ the hybridization function, such that A(cc;) = —iAo for a wide flat- 
band host [4] , with Ao = TrV^^phost ((^ = 0) the hybridization strength (a; = refers to the 
Fermi level, and V = V^k)- S(a;) = S^(a;) — iT,^{uj) denotes the interaction self-energy, 
which is merely defined by the Dyson equation implicit in equation (2.2). 

It is of course S(a;) that provides the conventional theoretical route to dynamics, 
via perturbation theory (PT) in U [4]. This approach is fine in principle, order by 
order in PT. But it is limited in practice, reflecting the inability of PT to handle strong 
correlations in general [4]. The practical difficulties arise because construction of Ti{uj) 
via conventional PT based on Hartrcc-modificd propagators, beginning with the static 
Hartree bubble diagram amounts to an expansion about the Hartree mean-field 
saddle point. And that single-determinantal saddle point is generally unstable to local 
moment condensation [3], reflected e.g. in the fact that the standard diagrammatic 
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resummations one might expect to be required to capture strong correlations — such 
as the sum of all particle-hole interactions in the transverse spin channel — give rise to 
well known divergences/non-analyticities in S(u;) [4,8,9]. 

In these circumstances the LMA simply recognises [8-10] that the natural mean- 
field saddle point about which to expand is unrestricted Hartree-Fock, corresponding to 
condensed local moments. Since H is invariant under a —a, this saddle point is now 
of course doubly degenerate, denoted hy a — Aor B corresponding respectively to local 
moments // = +\iJ,\ and — Accordingly, the full G{u;) is expressed as [8-10] 

G{uj)='^[GAa{io) + GBa{u^)] (2.3a) 

with self-energies flaai'^) (= + ^aaii^) with the static Hartree-Fock bubble). 
Note that equation (2.3a) is rotationally invariant, since the invariance of H under 
a <-> —a implies generally that GAai^) — Gs-ai^) hence that G{u!) is correctly 
independent of spin a. Direct comparison of equations (2.2) and (2.3b) also clearly 
implies a general relation (equation (3.4) of [9]) between the single self-energy S(a;) and 
the {Sao-(t^)}, enabling the former to be obtained directly from the latter. 

At this stage we emphasise the generality of the above considerations. The impurity 
propagator G{u), and hence spectrum D{ui), may be obtained either via the conventional 
single self-energy description embodied in equation (2.2), or via the two-self-energy 
(TSE) description inherent in equation (2.3): which one uses is a matter of choice, 
at least in principle. In practice, of course, the stability of the underlying mean-field 
(MF) saddle point arguably renders the TSE description a more natural choice. The 
standard diagrammatic resummations for the dynamical contribution EQ,c.(a;) to 'Eaa{'-^) 
no longer suffer from non-analyticities, and as detailed in [8-10] may be employed with 
impunity. This forms the practical basis of the LMA. Moreover the relation mentioned 
above between T,{uj) and {J^aai^)} may be employed to ascertain directly, and generally, 
the conditions on {^^^.(ci;)} under which ~ 0{u!^) as a; — >• 0, i.e. exhibits Fermi 

liquid behaviour. This is merely a matter of algebra, and as detailed in [9] the requisite 
condition is 

Sl(0) = EL(0) (= E^(0)) (2.4) 

referring exclusively to the Fermi level o; = 0. In an exact theory for the metallic AIM, 
this 'symmetry restoration' condition should be satisfied automatically (see also §3). 
For the LMA in practice [8-10] it is enforced self-consistently, and thereby determines 
the local moment magnitude (supplanting the usual gap equation for |//| at pure 
MF level) [8-10]. It is also central in being able to access the quantum phase transition 
to a degenerate local moment phase where such arises, as it does in the pseudogap 
AIM [12-14,16]. 

How well the LMA in practice captures single-particle dynamics can of course be 
tested by direct comparison to NRG calculations. For the metallic AIM, that has 
been considered in [10] (see also [14]). The resultant LMA scaling spectrum Di^uj) 



(2.36) 
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vs uj/uk in the strong coupling, Kondo regime was shown to be in excellent agreement 
with NRG results. However an issue evidently remains. Within the LMA the single- 
particle spectrum is expressed (equation (2.3)) as D{uj) = ^[Da^j^u) + Dbo-(u;)] (where 
Daa{^^) = —^^^Gacr{^))- It is this that has hitherto been compared to NRG results. 
The considerations above suggest however that there should exist an NRG counterpart 
of the spectral densities DAaii^) and DBa{^) inherent to a TSE description. It is this 
we now consider. The resultant Daa{^) can in turn be individually interrogated and 
compared to the detailed predictions for them arising from the LMA in practice. 

2.1. NRG spectra 

In Wilson's NRG method [1,2] the AIM is mapped onto a semi-infinite chain. This is 
diagonalised iteratively starting from the free impurity, and with a suitably truncated 
basis such that with increasing chain length essentially only the lowest lying states 
are renormalised. Spectral functions at each iteration N are calculated [5, 6] from 
the appropriate matrix elements connecting the NRG ground state and excited states, 
themselves related recursively to those of the previous iteration; and the spectrum 
for the whole frequency range is built up from the results for all iterations (roughly 
speaking, iteration determines the spectrum at frequencies lo ~ jJA.-i^-'^)/^ with D 
the conduction electron bandwidth and A the usual NRG discretisation parameter [1,2]). 
The resultant spectrum consists of a set of 5-functions with known weights at discrete 
frequencies, which is then broadened on a logarithmic scale to recover the continuum, 
see specifically equation (10) of [7]. This is the standard approach to calculating T — 
spectra via the NRG, and we follow it here. [An alternative method [7] focuses directly 
on obtaining the sclf-cncrgy S(c<j). Although not used here, we simply remark that it 
can also be extended to encompass the TSE spectral description considered below.] 

NRG spectra may be calculated from either the even set or the odd set of iterations 
(the fixed point, here exclusively strong coupling, is a fixed point of the square of the RG 
transformation [1,2]). The same spectrum D{uj) results in either shown explicitly 

below. There is however an important difference between even- A" and odd- A" iterations. 
In the former case, the NRG ground state is always a non-degenerate spin-singlet. For 
the odd- A" iterations by contrast we find that the NRG ground state (for any iteration) is 
a degenerate doublet, S — ^ (where we emphasise that 5' denotes the total spin angular 
momentum quantum number of the entire system). Denoting the Sz = +\ component 
of this doublet by 'A', and the = — | component by 'B', we can thus construct 
separate single-particle spectra (for cither spin a) from each of these degenerate ground 
states, denoted by DAa{^^) and DBaiy^) respectively; such that the total, normalised 
single-particle spectrum is given simply by D{uj) = ^[DAa{^) + -D_Bcr('-^)]- From the 
invariance of H under a <-> —a it follows that Dao-(<^) — Ds-ai^)', whence only (say) 
Dai{uj) and Dbi{uj) need be considered in general. And for the particle-hole symmetric 
AIM (ej = — -j) considered explicitly below, it follows further that Daa{^) = Da-a{—^^) 
{a — A or B), such that Db-\{uj) — Dai{uj) — Da-\{—u;) and D{u;) — D{—u;). 
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The essential point of the preceding discussion is clear. NRG spectra arising 
from odd-A^ iterations generate in effect a TSE description of single-particle dynamics, 
enabling both the D^ai^^) and the total spectrum D{u) to be determined. The even-A^ 
NRG spectrum by contrast, arising as it does from the non-degenerate NRG ground state 
in this case, corresponds to a conventional single-self-energy description of dynamics, 
from which D{uj) itself may again be obtained directly. 

3. Results: NRG and LMA 

To illustrate the above, figure [T] shows the resultant odd-A^ NRG spectrum Dai{uj) vs 
cu/Aq, for the symmetric AIM with a strong coupling interaction strength U = fZ/vrAo = 
12 (calculated with A = 2, retaining ~ 10,000 states per iteration). We consider first 
the insets, which show the corresponding total spectrum D^u) = ^[Dai{uj) + Db-\{uj)] 
from the odd- A iterations (solid line), as well as D{uj) obtained from the even-iV 
iterations (dotted line). The resultant symmetric spectra consist as expected of two 
high-energy Hubbard satellites, symmetrically disposed at cu ~ ±y, and the low-energy 
Kondo resonance. The first fact to note is that the D{u!) obtained from odd and even 
iterations are indistinguishable. This underscores the point made in §2: in an essentially 
exact approach, such as NRG, the two-self-energy and single-self-energy descriptions 
are fundamentally equivalent, and which one employs is a matter of choice. We also 
add (c/ §2) that the symmetry restoration condition equation (2.4), S^^(O) = E^^(O) 
(= for the symmetric AIM) is as expected satisfied within NRG, reflected in the 
fact that 7iAoDa^{uj = 0) = ttAqDai{uj = 0) = 1 (as required by the Friedel sum 
rule/Luttinger integral theorem [4], and satisfied in practice to < 1% accuracy in the 
present calculations). 

Turning now to the main panels in figure d the obvious feature is the intrinsic 
asymmetry in Da'\{uj). Only a lower Hubbard satellite is seen in Da-\{uj) (the upper 
satellite correspondingly arises in Dbi{uj) = D^|(— c<j)); consistent with the expectation 
that the A' state {S^ = +^) connects in the atomic limit (Vik = 0) to a purely |-spin 
occupied impurity, from which an |-spin electron may be removed but not added. The 
asymmetry in Dai{u}) is not moreover confined to the high-energy spectrum highlighted 
in the upper panel, which has an obvious counterpart in the atomic limit: the low-energy 
Kondo resonance (lower main panel) is likewise seen to be strongly asymmetric. 

To gain further insight into these results we now consider direct comparison with the 
LMA in practice, which affords an analytical handle on the single-particle dynamics. Full 
details of the underlying calculations are given in [8] for the symmetric AIM considered 
here (and in [9] for the asymmetric case). Here we focus purely on results arising. 

In figure El NRG and LMA results for nAoDA^iuj) are compared, vs uj/Aq (i.e. on 
an 'absolute' energy scale). The chosen f/ = 12 is simply representative of the large-?/, 
strong coupling behaviour of interest; similar results to those that follow have naturally 
been obtained for a range of interaction strengths. Let us first comment on the lower 
Hubbard satellite in DA^ii^j) (main panel), which clearly contains the vast majority of the 
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Figure 1. NRG spectrum ttAoDai{co) vs cj/Aq for U = 12. The upper panel shows 
-Da|(w) on 'all scales' including the single, high-energy Hubbard satellite and the 
Kondo resonance. The lower panel focuses on the asymmetric, low-energy Kondo 
resonance. Insets: the corresponding total spectrum D{lj) — i[I?^|(cj) + Z)b|(cj)] 
(= D{—u!)) obtained from the odd-iV iterations (solid line), and D{uj) obtained from 
the even-N iterations (dotted line); the two are in fact indistinguishable. 
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Figure 2. t:AqDa^{lo) vs lu /Aq for U = 12, obtained from NRG (solid line) and LMA 
(dotted line). The inset continues the comparison on a lower energy scale. NRG is 
known to overbroaden the high energy Hubbard satellite; this is reconstructed in figure 
6 below, where further comparison is made between NRG and LMA. 



spectral weight. The LMA sateUite is itself a Lorentzian, with a HWHM of 2Ao — twice 
that occurring in the pure mean-field limit (for the physical reasons explained in [8]), 
and which behaviour is believed to be asymptotically exact in strong coupling. With 
NRG by contrast the Hubbard satellite is well known to be overbroadened, due to the 
associated broadening on a logarithmic scale [5,6] (as occurs also, although ameliorated 
somewhat, with the method introduced in [7]). We will return to this below (see figure 
6), and show that the NRG results for the satellite are in fact consistent with the above 
Lorentzian behaviour. 

Aside from the high energy satellite issue, NRG and LMA results are seen from 
figure El to be in remarkably good agreement. That this persists down to much lower 
energy scales is shown further in figure El where Da^{uj) is shown in the top panel and 
the full symmetric spectrum D{u}) in the lower panel. Our aim now is to obtain a 
handle on the characteristic asymmetry evident in the Da^{uj) Kondo resonance. This 
may be obtained analytically from the LMA [10], formally for ^ 1 where ui = uj/ujk 
and the Kondo scale ujk is (here defined as) the HWHM of the Kondo resonance in 
D{ijj). On the positive frequency side, for = ^ 1, the LMA gives the asymptotic 
behaviour [10] 

TtAoD^tI^) = r4i , 1^12^1 (3-1) 

[-ln(a|aj|)]2 + 1 

with a ~ 0.7 a pure constant; while for negative frequencies by contrast, and —uj = 
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Figure 3. For U = 12, comparison between NRG (solid lines) and LMA (dotted 
lines) on the lowest energy scales appropriate to the Kondo resonance. Upper panel: 
7rAo£'AT('^) '^s f-^/Ao. Lower panel: the full spectrum ■k/S.qD{lo) vs oj/Aq, with 
D{uj) = i[i?^|(a;) + Z?BT(^)]- 



\uj\ > 1, 



ttAqD 



At 



[^ln(ap|)]2 + 25- 
Figure El compares NRG results for ttAo-DatI^) 



(3.2) 



with this predicted asymptotic 
behaviour. The agreement is seen to be very good, and in practice equations (3.1,2) 
already describe the spectra quite accurately for \u!\/u!k ^ 2 or so. We also add that 
the full forms of equations (3.1,2) are required for the agreement shown in figure EJ 
i.e. the behaviour in the |co'|/co'7^-range shown is not dominated by the ultimate high- 
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Figure 4. Comparison between NRG results for ttAoDat{u!) vs \uj\/ujk (solid lines) 
and the predicted asymptotic behaviour equations (3.1,2) (dotted lines). Lower curves: 
for positive frequencies, w > 0; upper curves: for uj < 0. 

frequency asymptotic behaviours of 7r^/[16 In^dcjl/cij/^ )] and 57rV[161n^(|a;|/a;K)]. For 
the full spectrum D(co') = ^[Dai{uj)+Db-\{uj)] = ^[Dai{uj)+Da-\{—uj)] the corresponding 
result is obviously just the weighted sum of equations (3.1,2) (with ultimate asymptotic 
behaviour of 37r^/[16 W(|c<j|/c<j/^)] that is exact for the s = | Kondo model [10]). It 
has been shown in [10] that this result for D{uj) itself is in excellent agreement with 
NRG results. We regard it as remarkable that the LMA provides an equally compelling 
description of the asymmetric Daa{^^)- 

The fact that the Daa{uj) are strongly asymmetric also means of course that the 
degenerate odd-iV NRG ground states 'A' and 'B' carry a local moment, equal and 
opposite for the two states, with magnitude given by |/i| = f^^du [D^|(u;) — Dai{uj)]. 
The resultant is shown in figure El as a function of U, and compared to its pure 
mean-field counterpart. The latter is given from solution of \fi\ = ~tan~^(^^) [15] 
(and in practice the LMA determined from symmetry restoration is very close to 
the MF value, exponentially so in strong coupling U ^ 1 [8]). Recognising that in 
strong coupling the Kondo resonance carries exponentially small weight and thus makes 
essentially no contribution to the resultant the leading asymptotic behaviour of 
the moment may be determined exactly from second order perturbation theory in the 
hybridization V^k- It is found thereby to be given by \^\ ~ 1 — A/ir'^U , which amusingly 
is precisely the leading result obtained from pure mean-field, and which as seen from 
figure El is indeed recovered correctly from the NRG calculations. 

We return finally to the issue of the overbroadened NRG Hubbard satellite, and 
show that the 'raw' NRG results are entirely consistent with a Lorentzian satellite 
of width 2Aq. To that end we calculate, as a sum of poleweights, the cumulative 
NRG spectrum FnrgI"^) — I-oo^^ Da-\{uj); the corresponding result for the pure 
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Figure 5. NRG local moment — J_^duj — _D^^((jj)] vs U — U/ttAq (solid 

line); compared to its mean-field counterpart (dotted line) which in strong coupling 
U ^ 1 gives ~ 1 — A/tt'^U. The latter result is in fact asymptotically exact, and is 
seen to be recovered by NRG. 

Lorentzian would be Fl = ^ ^^'^~^i '^2Ao ^ ~^ \ with ujq ~ — ^ the satelhte maximum. 
Writing FnrgI^^) = + 5F{u) with 5F{uj) thus defined, the NRG poles in Dai,{uj) 

contributing to 5F{ijj) are broadened in the usual logarithmic fashion, and added to the 
pure Lorentzian contribution (the net spectral weight below the Fermi level is of course 
preserved). If the NRG results are consistent with the Lorentzian, little deviation from 
this form should result. That this is indeed the case is shown in figure IHl and NRG 
and LMA results now agree well on all energy scales including the high-energy Hubbard 
satellite. 

4. Concluding remarks 

In this paper we have considered first a general issue: does a two-self-energy description 
of dynamics that is inherent to the local moment approach arise also within the 
numerical renormalization group method? The answer to that question is clearly 
yes — both it and a conventional single-self-energy framework arise naturally within 
NRG, according to whether odd or even RG iterations are considered. In consequence, 
explicit NRG results for the composite broken symmetry spectra underlying the two- 
self-energy description may be obtained: DAai,'^) and Db(t{^) such that D{uj) = 
^[DAa{^) + Dsai^^)] gives the total impurity spectrum (and with D{lj) coincident for 
both odd/even iterations). These in turn have been compared to results arising from the 
LMA as it is implemented in practice, and very good agreement found on essentially all 
characteristic energy scales from the high-energy Hubbard satellites to the low-energy 
Kondo resonance. 



Single-particle spectrum of the Anderson model 



11 



1.0 



0.8 



^ 0.6 

o 

t= 0.4 



0.2 




0.0 r~ 1 1 1 1 — — I 1 1 

-30 -25 -20 -15 -10 -5 5 

Figure 6. ttAoDa^{u)) vs lo/Ao for U = 12 obtained from the LMA (dotted line) and 
NRG (solid line), with the NRG Hubbard satellite obtained as described in the text. 
The agreement is now excellent on all energy scales. The original, overbroadened NRG 
spectrum is shown for comparison (dashed line). 



We also add that our essential conclusion is naturally not specific to the metaUic 
Anderson impurity model considered here: a two-self-energy description will arise for 
essentially any quantum impurity model, such as the pseudogap AIM [12-14, 16] (where 
the TSE description is in general a necessity and not a luxury), and more generally 
for lattice-based models such as the periodic Anderson model within the framework of 
dynamical mean-field theory [17]. 
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